Not to appear in Nonleamed J., 45. 



The dust properties of z ~ 3 MIPS-LBGs from photo-chemical models 

X.L. Fani'2 3.\ A. Pipino^ and F. Matteucci^-"^'^ 
ABSTRACT 

The stacked spectral energy distribution (SED) of Multiband Imaging Photometer 
for Spitzer (MIPS) 24^m detected Lyman break galaxies (MIPS-LBGs) is fitted by 
means of the spectro-photometric model GRASIL with an "educated" fitting approach 
which benefits from the results of chemical evolution models. The star formation rate- 
age-metallicity degeneracies of SED modelling are broken by using the star formation 
history and chemical enrichment history suggested by chemical models. The dust 
mass, dust abundance and chemical pattem of elements locked in the dust component 
are also directly provided by chemical models. Using our new "fitting" approach, we 
derive the total mass M,o,, stellar mass M*, gas mass Mg, dust mass Mj, age and star 
formation rate (SFR) of the stacked MIPS-LBG in a self-consistent way. Our estimate 
of M, = 8 X 10"' of the stacked MIPS-LBG agrees with other works based on UV- 
optical SED fitting. We suggest that the MIPS-LBGs at z ~ 3 are young (0.3-0.6 Gyr), 
massive (M,o, ~ IO^'Mq), dusty (M^ ~ lO^Mg), metal rich ( Z ~ Zq) progenitors of 
elliptical galaxies suff"ering a strong burst of star formation (SFR ~ 200Mo/yr ). Our 
estimate of = 7 x 10^ Mq of the stacked MIPS-LBG is about a factor of eight lower 
than the estimated value based on single temperature grey-body fitting, suggesting that 
self-consistent SED models are needed to estimate dust mass. By comparing with the 
Milky Way molecular cloud and dust properties, we suggest that denser and dustier 
environments and flatter dust size distribution are likely in high redshift massive star 
forming galaxies. These dust properties, as well as the different types of SFHs, can 
cause different SED shapes between high redshift star-forming ellipticals and local 
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star-burst templates. This discrepancy of SED shapes could in turn explain the non 
detection at submillimeter wavelengths, of IR luminous (L/« > IO^^Lq) MIPS-LBGs. 

Subject headings: galaxies: high redshift — galaxies: starburst — ISM: clouds — ISM: 
dust, extinction 



Introduction 



The last decade witnessed a tremendous increase in the knowledge of hig h redshi 
ing galaxies. Large samples, obtained with different selec tion techniques (see IShaplev 
review), yielded statistical information ab out masses (e.g 



ErbetaL 



2006b 



' t star form- 
201llfora 



Daddi et al. 



formation rates and dust attenuation (e.g. iHopkins & Beacomll2006l : iReddy et al.l 12008 
well as metallicities (e.g. lErb et al.ll2006a[) for z ^ 2 galaxies. 



2007 L star 



2010 ). as 



In a small subset of well studied objects, a more accurate determination of their chemical 
abundance pattern, d ust depletion, kinematics, was possible, d ue to lensing magnification (e.g. 



Mainline et al 



20111 : 



Genzel et al 



20091) and deep integral field observations (e.g. ILaw et al.l 1201 2l : lOnerucci et al 



201 Ih . 



as 



Taken together, these observables shaped our view of z = 2 - 3 redshift star forming galaxies 
relatively evolved systems, as witnessed by their metallicities (e.g. Erb et al. I l2006al: ICalura et al.l 



2009b . whose st ar formation rates are high due to the availability of large gas re servoirs as opposed 



to mergers (e.g. iGenzel et aLlBOlOl iRodighiero et al.l|201 It iKaviraj et al. 



2012). 



Among the various methods, the Lyman Break Technique (ISteidel & Hamiltonlll993h was 
particularly successful in yielding insight into high reds hift star forming galax ies. Lyman break 
galaxies (L BGs) are observed at diff erent redshifts: z~l (iBurgarella et al.ll2007h, z~ 1.4 -2.5 "BX" 
and "BM" Jsteidel et al.ll2004l) . z ~3 dsteidel et al.ll2003l) . z ~ 4,5,6 (IVanzella et al. II2OO9I) by means 
of the photometrical selection technique for rest-frame 912 A Lyman-continuum discontinuity in 
a desired redshift interval. To infer the properties of stars, interstellar medium (ISM) and dust, 
LBGs have been stu died in terms of stellar light, interste llar absorption lines and rest-frame UV 
continuum at z > 3 (ISteidel et all boOa Ivijh et al.l 120031) . One of the best studied LBGs is the 
gravitational lensed LBG MS 1512-cB58 (hereafter cB58) at redshift z = 2.7276, which seems 
to be a young system und ergoing rapid star formation and possibly destined to be an elliptical 



galaxy (IPettini et al.ll2002h . By means of cosmological hydrodynamical simulations. iDave et al. 
(|2000|) suggested that LBGs at z ~ 3 are the most massive galaxies at that epoch without rul 



ing out the merger-induc ed star-burst galaxies contribution. Using a chemodynamical model, 
Friaca & Terlevichl(ll999h sugg e sted th at L BGs at z ~ 3 a re the progenitors of present-day massive 
spheroids. iMatteucci & Pipinol (|2002|) and IPipino et al.l (|201 1|) suggested that LBGs at z ~ 3 can 
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be young small star forming elliptical galaxies by fitting the abundance pattern of cB58 by means 
of gal actic chemical evolutio n models. The stellar population properties of LBGs at z ~ 3 are stud - 



ied bvlSawicki & Yeel(ll998[) : IPapovich et all (1200 lh : lMagdis et all (120081 l2010ah : IPentericci et al 

hold . 



In the absence of accurate spectroscopic data available for large galaxy samples, the stellar 
properties, such as the age and stellar mass are estimated by spectral energy distribution (SED) 
fitting based on stellar population synthesis models where star formation history (SFH), metallicity 
and dust extinction are independent free parameters. Whilst some studies adopt a full treatrnent of 



radiative transfer (e .g. iTakagi et al.l 12003 



Siebenmorgen & Krugelll2007l : lRowan-Robinsonll2012 



Lo Faro et al.ll2013|) . in most cases empirical attenuation laws are usually adopted to deconvolve 
the eff'ect of the dust from the SED. This means that molecular cloud (MC) condition, and dust 
properties, such as dust mass, dust abundance and dust size distribution, are not directly studied . 
Only a few works attempt to fit the SEDs of LBGs from the UV to the IR (e.g. iMagdis et aPboiObl . 
201lb . However, separate models are adopted in these works to fit the UV-optical and the IR parts 
of the SED. Furthermore, the dust temperature and mass are usually estimated by a grey body 
fitting, and the IR luminosity are usually estimated with template SEDs from an observational 
library. Whilst these assumptions might be justified as a necessary to explore uncharted territory 
and certainly provide a good zero order estimate of the relevant physical quantities (masses, SFRs), 
we believe that more sophisticated models should be adopted to understand the physical processes 
acting in high redshift star forming galaxies at a finer degree of detail. 

To this purpose, in previous work (Matteucci & Pipino, 2001, Pipino et al., 201 1), w e focussed 

2000) and 
2009). We 



on the chemical abundance pattern of a well studied sing le LBG cB 58 (IPettini et a 



LBGs from two surveys: AMAZE (|Maiolino et al.ll2008h and LSD (|Mannucci et al 



self-consistently derived the SFH and the dust properties for these high redshift LBGs. Here we 
want to address in detail the dust properties of z ~ 3 LBGs. In particular, we will make use of 
the st acked SED of Multib and Imaging Photometer (MIPS) 24/zm detected LBGs (MIPS-LBGs) at 
z ~ 3 (|Magdis et al. koiObI) . since they allow for the first time an entire coverage from the UV to the 
radio wavelenghts. MIPS-LBGs, per se, are not a special class of objects. However they are likely 
the most massive, most dust rich, most rigorously star forming LBGs at redshi ft z ~ 3, following 



a SFR-st ellar mass relation with a similar slope to that found at redshift z=l (lElbaz et al.ll2007|) 
and z=2 (IDaddi et al.ll2008|) . but a higher normalization factor ([Magdis et al.ll2010a|) . The inferred 
properties of the stacked MIPS-LBG SED will represent the average massive LBGs at redshift z ~ 
3. The high SFR and infrared luminosity suggest that MIPS-LB Gs should be detected at submm 



wave l engths. However, this occ urs in a small fraction of cases dChapman et al.ll2000l : IShim et al. 



20071 : IChapman & CasevI 120091) . This non detection far-IR counterparts of MIPS-LBGs could 
be caused by higher dust temperature, different dust spati al distribution and lower SFR in LBGs 



than in submillimetere galaxies (SMGs) (see discussion in IChapman et al.ll2000l : lRigopoulou et al 
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20061 : iMagdis et al.ll2010b|) . The self-consistent SED model combined with chemical model will 
allow us to investigate this issue in more detail. 



In our approach (see details below), we study the 



the spectro-photometric model GRASIL (ISilva et al. 



for any galactic SED fitting (see reviews by iGawisei 



prope rties of stellar population and ISM by 



19981) . Parameter degener acy is a problem 



2009 



Walcher et al. 



2011). Therefore, our 



self-consistent chemical evolution model (jPipino et al.ll201 11) . which can reproduce most of the 
chemical properties of ellipticals of different mass, is adopted to reduce the degree of freedom 
as much as possible in our SED modelling by GR ASIL. The combinat ion of chemical evolution 
models and GRASIL has been already adopted by ISchurer et al.l (|2009|) to model SE Ds of differ 



ent mo rphological type galaxies using the chemical and dust evolution models from iCalura et al. 
(I2OO8I) . 



In this article, we aim at: 1) deriving the dust mass, and 2) investigating dust composition, 
temperature and grain size distribution of MIPS-LBGs. We will also assess the robustness of 
previously derived galaxy masses and SFR, compare MIPS-LBGs with other z ~ 3 LBGs, confirm 
the nature of LBGs at z ~ 3 from previous chemical models, and try to solve the challenge of non 
detection far-IR counterparts of MIPS-LBGs. The paper is organized as follows: in Sect. HI we 
describe our models and fitting approach; our results and discussions are presented in Sect. |3l and 
our conclusions are drawn in Sect. 51 Throughout the paper, we adopt a (0.7, 0.3, 0.7) cosmology. 



2. The models 



We combine the chemical evoluti on models for elli ptical galaxies (IPipino et al.ll201 11) and the 
spectro-photometric model GRASIL (jsilva et al.lll998l) to model t he SED of the stacked MIPS- 
LBG. The combination approach is fully described in lSchurer et al.l (120091) for galaxies of different 
morphological type. We direct the reader to above articles for equations and details. For the sake 
of convenience the two models will be briefly summarized in the next sections. In particular, the 
details of the our "fitting" approach will be described. 



2.1. The chemical evolution model 



The chemical abundance ratios versus metallicity are amongst the best indicators to constrain 
the SFH of a galaxy. However, there is no chemical data of th e average MIPS-L BGs at z ~ 3 
to constrain the SFH. We adopt the SFHs of elliptical models (IPipino et al.l 1201 1[) . which have 
successfully reproduced the chemical abundance properties of no rmal ellipticals (se e Pipino & 
Matteucci, 2004) as well as those of low mass LBGs at z ~ 3 (see IPipino et al.ll201 iL for detail). 
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In this way, the SFH, metallicity enrichment and dust evolution history of a galaxy are obtained 
from a self-consistent chemical evolution model. The impact of this assumption is reviewed in the 
discussion section. 

In these models, both the infall and the star formation timescale ar e assumed to decreas e with 
galactic mass, in order to reproduce the "chemical downsizing" (e.g. iThomas et al.ll2005|) . The 
initial conditions for ellipticals allow for the formation by either collapse of a gas cloud into the 
potential well of a dark matter halo or, more realistically, by the merging of several gas clouds. 
In any case, the timescale for both processes should be shorter than 0.5 Gyr, so that the ellipticals 
form very rapidly. The rapid gas assembly triggers an intense and rapid star-formation process that 
lasts until a galactic wind, powered by the thermal energy injected by stellar winds and SNe (la, II) 
explosions, occurs. At that time, the thermal energy is equal or larger than the binding energy of 
gas, and all the residual gas is assumed to be lost. After the wind, the star formation stops and the 
galaxies evolve passively. The evolution of the global metallicity as well as of 21 single chemical 
elements are studied in detail. Those elements are produced by single low and intermediate mass 
stars (AGB stars), SNe la, single massive stars (SNe II) and Type la SNe( white dwarfs in binary 
systems). The yields used in this paper are as follow s: 1) For single low and interrn e diate mass 
stars (0.8 < M/Mq < 8) we make use of the yields by Ivan den Hoek & GroenewegenI (Il997h as a 



function of metallicity; 2) For SNe la and SNe II w e adopt the empirica l yields by Francois et al. 



(12004^. These yields are a revised version of the IWoosley & Weaver (mi (for SNe II) and 
Iwamoto et al.l (119991) (for SNe la) calculations adjusted to best fit the chemical abundances in 
the Milky Way. We consider that the dust producers are AGB stars, SNe la and SNe II. We also 
take in to account the dust accretion and dust destruction processes in the ISM. Only the main 
refractory elements, C, O, Mg, Si, S, Ca, Fe, are depleted into dust, and we assume that stars can 
produce two different types of grains: i) silicate dust, composed of O, Mg, Si, S, Ca and Fe; and ii) 
carbon dust (graphite in the GRASIL language), composed of C. Since there is no AGN signature 
in MIPS-LBGs rest-frame UV spectra as well as in their SED, we do not include QSO dust in this 
paper. 



2.2. The spectro-photometric model 



Using the SFH, chemical and dust evolution of a galaxy extracted from the chemical evolu- 
tion model (e.g. as shown in Fig[T]and F'ig. |2(a)(|2(c)|), we synthesize the SED of that galaxy with 
the spectro-photometric code GRASIL (^ Silva et aL 1998 ). Briefly, emissions form stell ar popula- 
tions are calculated by evolutionary population synthesis technique feressan et al. 1994 ) us ing th e 
chemical evolution model and the Padova simple stellar population model (IBertelli et al.lll994|) . 
Younger stellar generations are more aff"ected by dust obscuration in their birth place (molecular 
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Table 1 : Chemical model parameters and results. 



model 
name 


(Mo) 


(kpc) 


V 

(Gyr-') 


T 

(Gyr) 


tgw 

(Gyr) 


SW 
yields 


SNe la and SN 11 
yields 


dust elements 


M310 


3 X 10'" 


2 


3 


0.5 


0.8 


V&G 


Francois 


C, Si, Fe, Mg, O ,(S, Ca) 


Mil 


10" 


3 


10 


0.4 


0.7 


V&G 


Francois 


C, Si, Fe, Mg, O, (S, Ca) 


M511 


5 X 10" 


6 


15 


0.3 


0.7 


V&G 


Francois 


C, Si, Fe, Mg, O, (S, Ca) 



M;„„: the total mass of the galaxy at galactic wind time tg„. 
Reff- the effective radius of the galaxy, 
v: the star formation efficiency. 
t: the infall time scale. 

fgH- '■ the galactic wi nd time. It is determined only by M/„m (and hence v and t). Note that it i s shorter than i n 
IPipino et alj |201 lb . This is because the (g^. refers to the central region of the galaxy in Pi pino et al 1 bo 111) , 
while it r efers to the whole galaxy in this work. This is a feature of the outside-in formation formation 

scenarios jPipino & Matteuc"dll2004h . 

Single low and intermediate mass stars yields (SW yields) are from lvan den Hoek & GroenewegenI l ll997h 

(V&G). 

SNe la and SNe II yields are from lFrancois et al] j2004l) . 

C, Si, Fe, Mg, O are the dust elements adopted in GRASIL, while S, Ca are also the dust elements in 
chemical models. 




Fig. 1. — The SFR versus time in different mass galaxies. M511, Mil and M310 represent 
5 X 10'\ 10^' and 3 x 10'*' Mq galaxies, respectively. The red horizontal line corresponds to a 
SFR=2OOM0/yr, which is the estimated value for the average MIPS-LBGs. 
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(a) Total (red solid line), graphite (black dot (b) Dust to gas ratio of lO^M© galaxy as a func- 
line) and silicate (blue dashed line) dust mass of tion of galactic age. 
IO'^Mq galaxy as a function of galactic age. 



0.030 



N 0.015 



o.ooc 




galactic age [Gyr] 

(c) Metallicity of 10"Mq galaxy as a function of 
galactic age. 



Fig. 2. — Output of the chemical evolution models adopted by GRASIL. Values in the star forming 
phase are shown, since we only focus on this phase. 
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clouds) than older ones in diffuse ISM. The effect of the age selective extinction of stellar pop- 
ulations, f or the first time, was modeled by a parametric approach adopting the parameter tesc in 
GRASIL (|silva et alJll998|) . This is described in GRASIL assuming that the fraction of starlight 
radiated inside the clouds by stars is a function of the star age. In practice, if tesc is the timescale 
for the process, 100% of the stars younger than tesc are considered to radiate inside the MCs, and 
this percentage goes linearly to 0% in Itesc- The dust composition in GRASIL consists of graphite 
and silicate grains, and Polycyclic Aromatic Hydrocarbons (R^H) molecules, with a distribution 
of gr ain sizes for each comp osition. The optical properties of silicate and graphite gra ins are taken 
from lLaor & Drainel (|l993h . The optical properties of RA.H molecules are taken from lDraine & Li 



(|2007|) . Dust abundance and dust composition relative abundances are obtained directly from the 
chemical evolution model. For the molecular clouds a full radiative transfer calculation is per- 
formed. The radiative transfer of starlight through dust is computed along the required line of 
sight yi elding the emerging SEP dGranato & DaneselfljgS) based on lambda-iteration numerical 
method ( Efstathiou & Rowan-Robinsonlll99dlCollison & Fix|[r99 lL However for the diffuse dust, 
the effects of scattering are only approximated by assuming an effective optical depth related to 
the true absorption and scattering optical depths by: Teff = \jabs{Tabs + Tsca)V'^ and assuming that 
there is no re-absorption of the radiation emitted by the cirrus. 



2.3. The difference between the chemical and spectro-photometric model 



The adopted chemical evolution models contain only one gas phase whereas GRASIL takes 
into account two gas phases: the molecular clouds and diffuse gas. We therefore assume that 
the chemical and dust-to-gas ratio and dust abundances are the same for those two gas phases. 
Ca and S are not contained in GRASIL's dust composition. However, those two elements give 
a negligible contribution to dust-to-gas ratio. In this paper, we will only consider the following 
refractory elements: C, Si, Fe, Mg and O. That is, quantities as the total dust mass will be computed 
neglecting other elements (see |2(b)| ). The effect caused by the difference between ISM (ISM= 
gas -I- diist) me tallicity and gas phase metallicity on SED is also negligible. The Salpeter IMF 
(|Salpetenll955h is adopted in the two models. The baryonic matter (i.e. stars plus gas) is assumed 
to follow the Ijaffd (|1983|) spatial distribution in the chemical evolution model, while GRASIL 
adopts the King profile ( King|[l966 ) spatial distribution. The spatial distribution, which determines 
the potential well, only directly affects the galactic wind time in chemical evolution models (see 
dis cussion in section 12. II). The geometry effect is not important in the elliptical galaxy SED (e.g. 



see 



Silva et al.lll998l : |Piovan et al.ll2006bl . for detail). The effects of various parameters on the SED 



will be discussed in further work in detail. In this paper we only describe the general effects of 
parameters on deriving physical properties of galaxies. 
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Table 2: Results from SED fitting for the average MIPS-LBG 



Models 


MW" 


MIPS ^ 


other works 


IR luminosity: log( L/r/Lq) 


12.13 


12.11 


12.2r 








12.44^ 


Stellar mass: M* (Mq) 


8 X 10'° 


8 X 10'° 


7.9 X 10'° " 


age (Gyr) 


0.5 


0.5 




SED-derived SFR (Mo/yr) 


200 


200 




IR-derived SFR 


233 


222 


275« 


Dust mass: (Mq) 


7 X 10^ 


7 X 10^ 


5.5 X 10^'' 


Observational frame 850 fxm flux density : f^5Q{mJy) 


0.16 


0.34 


1.36' 


UV Slope: /3 


-1.6 


-2.1 




Extinction: E(B - V) 


0.213 


0.127 


0.1 ^ 


Luminosity weighted metallicity: < Z >/„,„ 


0.017 


0.017 




ISM metallicity: Z 


0.02 


0.02 


0.02*^ 


Dust to gas ratio: Dg 


0.0025 


0.0025 





a: 10"Mo model (Mil) at 0.5 Gyr with WM-like parameters (see Table H and |4l) 
b: 10"Mo model (Ml 1) at 0.5 Gyr with MIPS paramete rs (see Table |3] and g]) 



c, g and i fro mlMaedis et al. 



(2010bl 



d and h from lRigopoulou et al.l (I201(]t) 
e, f and k from Magdis et al. ( 2010a|) 



j from Magdis's private communication 
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Table 3: Adopted values of MC in GRASIL 





MW 


MIPS 


t 1 

'■esc 


2 


200 


Jmc 


0.5 


0.5 




10^ 


10^ 


n 4 


40 


16 



^ Timescale for the evaporation of MCs, in Myrs. 
^ Fraction of gas content in the MCs. 
^ Total gas mass in each MC, in Mq. 
Radius of each MC, in pc. The value of R^c in the MW is just to respond to the lower 
dust-to-gas ratio in the MW than in the star-burst galaxies (see text in Sect. 13.21) 



Table 4: Parameters for size distribution of the dust cor nponents (Eg. [H a nd O in the MW and 
MIPS-LBGs. The MW size distribution is that derived by lSilva et al.l (Il998|) to match observations 
from the MW. The MIPS-LBGs size distributions are calculated in this paper to match SED. 





MW 


MIPS LEG 


PAH 
X 


3.3 X 10-25 cm"/// 


5 X IQ-^^cm^-^IH 


Graphite 






^min (^) 


8 




at (A) 


50 


8000 


O-max (A) 


2500 


22500 




-3.5 


-3.5 




-4.0 




Silicate 






ajnin (A) 






Ub (A) 


50 


800 


O-max (A) 


2500 


12500 




-3.5 


-3.5 
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2.4. The fitting approach 

In this work, wit h the help of the SF H, chemical pattern and dust properties from the chem- 



ical evolutio n models (iPipino et al.ll2011l) and pan-spectral energy distribution data of the stacked 
MIPS-LBG (iMagdis et alj 1201 Obi) , we break the SFR-age-metallicity degeneracies for SED mod- 
eling. Basically, the flux derived by median stacking analysis represents the average properties of 
many undetected individual objects. We will not give "ad hoc" parameters but the general proper- 
ties of the average MIPS-LBG by modelling the stacked SED without taking into account the filter 
widths in SED fitting. We will adopt an "educated fitting" approach to estimate the physical pa- 
rameter of the galaxy. This approach is guided by answering these questions: 1) Which parameter 
dominates the overall-level SED of a galaxy? In other words, from a galaxy SED, which estimated 
property is most reliable? 2) How do other parameters afl'ect the SED in detail, such as the shape 
and the peak of some parts of a SED? By means of this approach, we use SFHs suggested by self- 
consistent chemical evolution models. Moreover, chemical and dust properties, such as metallicity 
and dust mass, which are the basic ingredients in SED modeling, are given by the same chemical 
evolution model which gives the SFH. The "fitting" ap proach is the followin g: first, since the mass 



is the most robust parameter in SED fitting (e.g. see iShapley et al.ll2005|) . the total mass of the 
average MIPS-LBG is estimated by comparing the predicted overall-level of SED with data, and 
confirmed by further fitting steps. This treatment effectively reduces the computation time, since 
we do not need to test the MC and dust parameters in GRASIL to estimate the total mass ( Step 
One, see FigO. Once the galaxy is selected by total mass estimate, the SFH, chemical and dust 
evolution history are given by the chemical evolution model with that total mass. Second, the age 
of the average MIPS-LBG is estimated by comparing the total mass-selected overall-level of SED 
with data. At this step, the SFR, metallicity, stellar mass, dust mass of the average MIPS-LBG are 
derived by the SFH adopted in the chemical evolution model ( Step Two, see FigU). Third, the 
"best-fitting" dust parameters (shown in Table H] and ID are derived by min-;^f^ method based on 
a grid of testing parameters in GRASIL (see an example of grid models in Fig. [S]). At this step, 
the galactic parameters estimated in previous steps are confirmed by the "best-fitting" SED ( Step 
Three, see Fig|6l). 

To fit the SED, we have tested in GRASIL many SFHs origi nating from ellipti cal models. 
These models cover a large mass range (from lO^M© to IO^^Mq, see Pipino et al. 201 l l for detail). 



The parameters and results of three examples of tested models with total mass 3 x IO^^Mq (M310), 
lO^^Mo (Mil) and 5 x IO^Mq (M51 1) are shown in Tabled SFHs of these three models and the 
estimated SFR of the average MIPS-LBG are shown in Fig. [TJ Dust mass, dust-to-gas ratio and 
metallicity as a function of galactic age of a typical elliptical of IO^^Mq (the "best-fitting" galaxy, 
see below) adopted in the GRASIL are shown in Fig. |2(a)P(c)[ respectively. 



The free parameters in GRASIL, which will be investigated by min-;^^^ method in this work. 
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are the ones related to MCs and dust size distribution. There are four parameters (in Table |3) 
related to MCs: 1) the time scale t^sc of young stars escaping from their birth places (MCs), 2) 
the fraction of gas content in MCs f^c, 3) the mass of a single MC M^c, 4) the radius of a single 
MC Rmc- Note that, the predicted SED depends on M^c and R^c only through the combination 
Mmc/Rlic, which is the true free parameter. The dust size distribution for each dust component is 
described by the power laws (Eq. [Hand Eq. |2]in Sections 3.2 and 3.3). The parameters of graphite 
and silicate (in Table Sj) are the size lower limit (a„„„), size upper limit (a,„ax), the connection point 
of the two power laws (a/,) and the ones related to slopes (fii and /32 in Eq. O. The parameter 
X of PAH size distribution controls the abundance of PAH. The normalization parameters, which 
control the dust abundance and dust component relative abundances are given by th e chemical 



evolut ion model (see Sect. 13.21 and [33] for detail). Other parameters are the same as in [Silva et al. 



([1998 1. The standard (Milky Way, MW for short) and our "best- fitting" model parameters adopted 
in SED modeling are shown in Tables [2] and [H Effects caused by various IMFs, simple stellar 
populations and dust optical properties will not be discussed in this paper. 



3. Results and discussion 

The zero order estimates of the relevant physical quantities (masses, SFRs) are typically de- 
rived by global or integral properties, such as overall-level SED and IR luminosity, therefore es- 
timates of those quantities are robust. Other quantities (such as UV slope and the predicted flux 
at a particular band) depend on the shape of SED, and could be different. Table [2] shows the main 
fitting properties of the average MIPS-LBG by different models. These models fit the overall-level 
SED, and predict consistent IR luminosities. Stellar mass and IR-derived SFR in all models are 
comparable, while other quantities are not necessarily consistent with each other. We will discuss 
these properties in more detail below. 



3.1. Stellar populations 



The overall-level rest-frame SED, like the value of the luminosity density or flux density, 
of a galaxy is dominated by its total galaxy mass. Particularly, the total stellar mass dominate 
the UV-optical part of the SED and total dust mass dominate the IR part of the SED. In general, 
stars, such as AGB stars and SNe, are the main sources of dust in a galaxy. A more massive 
star-formin g galaxy contains m ore stars as well as more dust (e.g. see the case of young elliptical 
galaxies in [Pipino et al.[[201 1|) . The more massive star forming galaxy will show higher overall- 
level rest-frame SED. Since our estimations are based on the UV-radio overall-level SED fitting, 
our results of the stellar population are derived by a complementary approach comparing to the 
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noraial approach, which is usually based on optical to near-IR SED fitting. 

We compare the predicted SEDs of different massive galaxies at different ages (see an example 
of different massive galaxies at moderate age 0.5 Gyr in Fig. |J| Since short and intense star 
formation histories are yielded by chemical evolution models, the parameter space of age is not 
large (age < O.lGyr). From Fig. |3]and Fig. |4l we can estimate that the galaxy with total mass 
~ IO^^Mq could reproduce the overall-level SED of the average MIPS-LBG. Note that, this mass 
is the total mass of an elliptical galaxy at the wind time given by the chemical evolution model. In 
fact, after the wind, no more star are formed and the galaxy evolves passively. 

From the modeling point of view, the total stellar mass is dominated by SFH, IMF and the age 
of stellar populations. Age is the only free parameter of stellar populations when SFH, chemical 
enrichment history and IMF of a galaxy are fixed. With a given SFH in a chemical model, the SFR 
is a predictable parameter. By comparing with the value eatimated from observa tional indicators . 



SFR could constrain the age of a galaxy. UV and IR indicators for SFR (e.g. lKennicuttill998[) 
have been adopted in the literature for high redshift galaxies. Besides, depending on stellar popu- 
lation synthesis models, the UV-derived SFR depends on the treatment of extinction for deriving 
the extinction-correct fluxes, and the IR-derived SFR depends on the SED library for converting 
observed flux (e.g. the flux at 24 fxm) to total (8 - lOOOyum ) IR luminosity L[r. This approach 
has the drawback of assuming the low-redshift template SEDs accurately represent the SED of 
high-redshift galaxies. For example, the excessive high 850 /im fluxes predicted by an empirical 
UV-based relationship and the unreasonable high L/« estimated by MIR flux-IR luminosity cor- 
relation in SED t emplates have caused the p roblem of non-det ection of z ~ 3 LBGs in 850 fxm 



band, reported in IChapman et al.l (|2000h and IShim et al.l (|2007l) respectively. Additionally, vari 



ous dust intrinsic properties strongly affect the UV-FIR fluxes (see Fig. [5]). Therefore, we need 
other indicators to estimate the SFR in our SED fitting approach. In GRASIL the non-thermal 
radio emission, which is not affected by dust intrinsic properties, is proportional to the SN II rate. 
Since the life of SNe II is very short (< 30 Myr), the SN II rate is proportional to the SFR. So 
the rest-frame radio flux should be a good estimator of SFR. This is reflected by the fact that both 
the radio flux and SFR decrease (shown in Fig. [T]and|4l) at the age from 0.2 to 0.6 Gyr. For the 
lO'^Mg model, the phase with SFR ~ lOOM^/yr, corresponding to the age of 0.3 to 0.6 Gyr (see 
Fi g[D, fits the rest-frarn e radio flux well, and this SFR agrees with the estimation (~ 250Mo/yr) 



of iMagdis et al.l (|2010bh . Since we have already selected the 10^ 'Mq model based on overall-level 
SED, the SFR derived here benefits from the merits of radio-derived and SED-derived approaches, 
and is remarkably consistent with the IR-derived SFR (see Table [2] and discussion below). 

In our approach, the age is estimated by fitting the overall-level SED and the radio flux, 
namely it is a luminosity weighted age confirmed by "observed" SFR. The SFR and total mass 
suggest that the ages of MIPS-LBGs should be ~ 0.3 - 0.6 Gyr (see Fig. H), and younger than 
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0.7 Gyr since the LBGs are star-forming galaxies. Note that, the estimated SFR = IQQM^Iyr is a 
typical SFR in the age range of 0.3 to 0.6 Gyr. The SFR is lllMQiyr and UQM^Iyr at 0.3 and 0.6 
Gyr for the IO^Mq model, respectively. The total stellar mass is ~ 6 x 10'°Mq and ~ 1 x IO^'Mq 
at 0.3 and 0.6 Gyr for the 10' ^Mq model, respectively. We adopt 0.5 Gyr for the next fitting step. 
This age corresponds to SFR ~ 18 5Mo/yr and total ste llar mass ~ 8 x 10'"A/o, which agrees 
with the value (~ 7.9 x W^Mq) of Magdis et all (llOlObh . Our age estimate 0.3 - 0.6Gyr is not 
consistent with the median age ~ \Gyr derived by synthesis population model in iMagdis et al. 



(|2010a|) . In particular, our elliptical model does not allow a LBG to be older than \Gyr. It is 
however worth reminding, that while the star formation history predicted in chemical evolution 
models is modulated by the gas infall and self-consistently quenched due to SN-driven wind, this 
does not happen in the template star formation histories adopted in stardard SED fitting, and this 
might cause the disagreement in the age estimate. 



3.2. Dust properties 



As molecular clouds (MCs) are the birth place of stars, studying properties of M Cs at high 



redsh ift is important for our understanding of the star forming in galaxies (see review by iRiechers 
201 1|) . In this paper, we assume that all the MCs have the same mass M,„c and spherical radius 



Rmc- MCs make much more contributions in extinction than diffuse ISM in a star-burst galaxy with 
the fractionary gas content in the MCs /^c ^ 0.5. The parameter tesc shown in Table |3] controls 
how long the young stars remain in their birth clouds and roughly controls the fraction of extincted 
stellar light, therefore it controls the total level of rest-frame UV-optical SED and the slope of the 
rest-frame UV-optical SED. The optical depth of dust in a MC is determined by the dust-to-gas 
ratio 5 times M,„c/Rmc- This parameter moderately controls the total level of rest-frame UV-optical 
SED and also aff"ects the slope of rest-frame UV-optical SED. The "best-fitting" parameters of 
MCs for the average MIPS-LBGs are shown in Table |4l The larger value of the tg^c and the smaller 
value of the R^c compared with the MW one^ imply that the MCs in the average MIPS-LBG are 
likely to be in mor e dense dusty environments. The sartie trend is found in loca l and high redshif t 
starburst galaxies dsilva et al.lll998l: Ischurer et al]l2009[ Iswinbank et aPbouk lYan et all (I2OI0I) 
suggested that the average molecular gas mass of 24//m detected ultra-luminous infrared galaxies 
(ULIRGs, Ur > lO'^Lo) at z = 1.6 - 2.5 is 1.7 x 10'°Mq by observed CO J=2 ^ 1 or J=3 ^ 2 



^The values of MC parameters of the MW are the same as in spiral models in lSchurer et al.l (120091) . except for R^c- 
The dust-to-gas ratio is given by the chemical evolution models and not modified by GRASIL in this work. Since 
we do not use the spiral SED models, we adopt a larger value of Kmc for the MW responding to the effect of lower 
optical depth caused by lower dust-to-gas in the MW than in star-burst galaxies. In any case, the t^sc is the dominant 
parameter. 
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Fig. 3.— Rest-frame SEDs of different models. M511, Mil and M310 represent 5 x 10", 10" 
and 3 x lO'^'Mg galaxies with star burst SFH at 0. 5 Gyr with the MW d ust parameters shown 
in Table |3] and Table HI respectively. Data are from iMagdis et al.l (|2010bl) : for the SED we use 
the median UGR+BViJK+IRAC+MIPS24 photometry of MIPS-LBGs, and the values (or upper 
limits) derived from stacking 100 /im, 160 fxm AzTEC and radio, divided by 1 + z. Downward 
triangles are up-limits. The MIPS data is indicated with a red ponit. At the first fitting step, the 
10" Mq model is empirically selected for next fitting step. 




Fig. 4. — Rest-frame SED of a 10"A/o galaxy at different ages with the MW dust parameters as 
in Tableland Table 51 The parameter space of age is < 0.7 Gyr. The model with age 0.3 - 0.6 
Gyr could fit the overall-level SED. At the second fitting step, the lO^M© galaxy at 0.5 Gyr is the 
model adopted for further fitting step.wo This model predicts compara ble SFR and s t ellar rn ass 
with other people's work (see text for detail). MIPS-LBGs data are from lMagdis et al.l (|2010b|) as 
described in Fig |3] 
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emission. Our total gas mass Mg in lO^^M© galaxy at 0.5 Gyr is 1.6 x IO^^Mq. With fractionary 
gas content in the MCs /mc = 0.5, the total mass in molecular clouds in our "best-fitting" model is 
0.8 X IO^^Mq. However, the lack of emission lines makes it difficult to break the degeneracy of the 
MC conditions and dust size distribution. 

Once the emission from stellar populations is set, the peak in the emission in far-IR (FIR) is 
dominated by the total cold dust mass. The total dust mass A/j is ~ 7 x IO^Mq in "best-fitting" 
1 0^^ Mp, model at 0.5 Gyr It is about a factor of eight less than the value ( M4 = 5.5 ± 1-6 x lO^Mp 



in 



Rigopoulou et al. 20ld ) based on single temperature grey-body fitting ( Hildebrand 1983 ). This 



may partly be caused by three facts: 1) the uncertainty in the galaxy mass estimate, therefore the 
dust mass value in this work. Since the dust mass is a consequence of the galaxy evolution in 
our chemical models, the estimated total galactic mass aff"ects the total dust mass. For example, 
a 2 X 10"Mo galaxy § at 0.5 Gyr produces a dust mass ~ 2 x IO^Mq. 2) The well known 



degen eracy of dust temperature and slope fid in single temperature grey -body fitting ([Blain et al 



20031). A hi gher Pd will resul t in lower dust temperature derived form single temperature grey- 
body fitting (|Sajina et al.ll2006|) . therefore a higher dust mass will be estimated. 3) The uncertainty 
of single temperature grey-body fitting parameter, suc h as the rest-frame dust mass absorption 
coefficient k (e.g. a factor ~ 7 at 800 ju m estimated by ^ Hughes et al. 1997 ). Note that, the single 
temperature grey-bod y fitting could predict a higher IR flux than the value by full SED fitting. An 
example is shown in IMagdis et al.l (l2010b|) . The value of the far-IR emission peak predicted by 
the single temperature grey-body fitting shown in their Fig. 4 is larger than 1 mJy, while the peak 
value shown in their Fig. 2 is smaller than 1 mJy predicted by SED fitting. 

Prominent PAH features are observed in MIPS-LBGs ( Huang et al.ll2007 ) . In star-burst galax- 
ies, the mid-IR (MIR) flux is contributed by small hot dust and PAH (see e.g. iLaurent et al.ll2000l 
and references therein). The abundance of PAHs of our model is calculated from the chemical 
composition of the dust as predicted by the chemical evolution model, and their abundance is pro- 
portional to the total abundance of carbon molecules in the dusty component of the ISM. The PAH 
size distribution adopted in GRASIL is : 



^" V -3 5 -1 

— = Xa cm . 
da 



(1) 



PAH is needed in the "best-fitting" model to fit the average MIPS-LBG MIR flux (see FigUl). 
However, we cannot put more constraints on the PAH without emission l i nes. T he "best-fitting" 
model gives X = 5x \0~^^cm^-^/H (X = 3.3 x IQ-^^cm^-^/H in lsilva et all (1 19981) '). 



deviation of estimated total galaxy mass from the "best-fitting" value within a factor of two is normal and accept- 
able in all SED fitting work. 
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3.3. Dust intrinsic properties 



For many years, many efforts have been spent on investigating the du st size distribution 



through three approaches: i) theoretical approach (e.g. 



tion curve approa ch (e.g. iMathis et al.l 119771 IWeingartner & Draine 



Bimstiel et al. 



Zubko et al.ll2004l) andiii) by means of the SED approach (e.g. 



2004 



201 



2001 



Carciofi et al.l 



I, ii) fitting extinc 



Clayton et al. 



20041: 



2003 



Takeuchi & Ishii 



Piovan et al.ll2006ah . However, the nature of dust size distribution is still not clear, especially 



for in high redshift objects. 

The position of the peak in FIR is mainly affected by the total dust mass, dust component 
relative abiindanc es and dust size distribution, but not by the stellar populations (see Figj3] and 
Slater et al.ll2011|) . With the help of chemical evolution models, which predict the total dust mass 
and dust component relative abundance, we can study the dust size distribution in a more reliable 
way. The graphite and silicate grains dust size distribution strongly affect the slope of the rest- 
frame UV-optical SED. For the sake of simplicity, a simple power low dust size distribution form 
is adopted in this paper. The dust size distribution for the / (graphite or silicate) composition 
follows a broken power law defined with a threshold ab- 



dtii 
da 



ifah < a < Umax, 
'^min < a < Ub- 



(2) 
(2') 



For the sake of simplicity, dust in MIPS-LBGs only use Eq. |2](see Table|4]). The dust mass m, and 
dust-to-gas ratio 6 are provided by the chemical evolution model: 



m,- 



6 = 



4na^Pi drii 
3nH da 

niH 



da. 



(3) 



(4) 



where hh, mn and p, are the gas number density, gas mass and grain mass density, respectively. The 
normalized parameter A, is calculated by Eq. [2l [31 IH Comparing with the dust size distribution 
in the MW (see Table |4]), a larger amount of dust with larger size is needed to fit the slope of 



Tanaka et al 



(e.g. iRodmann et al 



rest-fram e UV-optical SEP (see Fig. [5] and 
size (e.g. 



l.ll2005L Bimsti eletal. 
20061 : iRicci etallbon 



Fig. [6]). With the theoretical indication of a few /urn 



2009 b and observational suggestion of up to cm-sizes 



20101) . our estimation of maximum dust size amaA^=2.2 fim 



is acceptable. 

With indications from both the peak and the position of FIR and the slope of rest- frame UV- 
optical SED, we suggest that the dust size distribution of the average MIPS-LBG is top-heavy 
(flatter) comparing with that in the MW. The reliability of this trend for all high redshift LBGs 
needs more data to provide more reliable constraints for the models in the future. 
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3.4. The impact of the adopted SFH 

The adopted SFH is the main input ingredient which produces different predictions on de- 
tailed physical processes by means of population synthesis modeling. By the standard SED fitting, 
we cannot rule out other SFHs than the ones adopted. For example, Fig. |7] shows two mod- 
els adopting a self-consistent star formation history (SFH of Mil case in FiglH SSF for short) 
and a constant star formation history (CSF). These two models can produce similar overall-level 
SEDs with fine-tuning of the ages corresponding to comparable stellar masses. In this paper, to 
improve the reliability of the adopted SFH for the z~ 3 MIPS-LBGs, we adopted chemical model 
constrained SFHs, which have been successfully reproduced the chemical data of several z~ 3 
LBGs. 

The stellar mass derived by Elliptical-type star-burst SFH in this work is consistent with other 
works, although we used a difi'erent SED model and fitting approach to fit the SED. Again, this is 
because the stellar mass is a more robust p arameter than other parameters, such as age, which have 
been fo und in oth er SED fitting works ( e.g. Papovich et al.lboOl : Shapley et aklllOOS : Magdis et al. 



2010a ^. Recently, Lo Faro et al. ( 2013 1) found discrepancies between their results and those based 



on optical-only SED fitting for the same objects. By fitting observed SEDs with their physical 
model, they found higher stellar masses (by AM* ~ 1.4 and 2.5 dex) resulting from higher extinc- 
tion (the average total extinction in the rest-frame V-band Ay ~ 0.81 and 1.14 ) for z ~ 1 and z ~ 2 
dusty normal star-forming (e.g. SFR< lOMojr"') galaxies, respectively. With only one average 
z ~ 3 dusty star- forming galaxies in our paper, it is hard to test whether there is also a same system- 
atic discrepancy, and the discrepancy is mainly caused by extinction for z ~ 3 dusty star-burst (e.g. 
SFR> lOOMojr"') galaxies or not. We would like point out that the dusty star- forming galaxies 
at difi'erent redshifts are not necessarily linked from an evolutionary point of view, therefore the 
types of SFH and the evolution phases in which galaxies are observed might be difi'erent. There- 
fore, one might get comparable stellar mass by the full SED fitting and optical-only SED fitting 
with different ages and extinction values (laws). For example, by playing with dust parameters in 
the CSF case shown in Fig. U\ we would expect to fit the optical-only SED at 2 Gyr estimating 
a stellar mass ~ 9 x 10^° which is comparable with the stellar mass ~ 8 x 10^° at 0.5 Gyr in the 
"best-fitting" model by fuU-SED fitting. 

Unlike the stellar mass, our estimated age (~ 0.5 Gyr) is younger than other works (~ 1 Gyr). 
Our intense star-burst SFH allows the galaxy to form enough stars and dust to produce the overall- 
level SED at that age. Dust makes the key link between the UV-optical SED and IR SED. In the 
standard SED fitting approach, dust mass cannot be derived, therefor it does not depend on the 
adopted SFH. In our approach, dust mass is a consequence of galaxy and dust evolution in the 
chemical models and depends on the SFH. By assuming reasonable detailed dust properties, the 
dust in our model interacting with stellar light, reproduces the SED of the average MIPS-LBG. The 
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value of the dust mass in our work is about a factor of eight less than the value in lRigopoulou et al 



(120 10|) based on single temperature grey-body fitting (see discussion in Sect. I3.2I ). Note that, our 
predictions on dust environment and dust intrinsic properties correspond to "best-fitting" SED of 
the average MIPS-LBG, and further work/data about these MIPS-LBGs are needed to confirm this 
result. 

Clearly, the high SFR shown in the LBGs could be triggered by galaxy-merger, which means 
complectly different SFH and physical processes. This scenario of galaxy evolution is not con- 
sidered in this work. We hope that more data in the future, such as the galaxy morphology and 
chemical abundances, will reveal the nature of all LBGs. 

Having discussed the limitations in our approach, we will compare our "best-fitting" model for 
the average MIPS-LBG with other objects which show similar observed features in some aspects. 



3.5. Comparison with lower mass z~3 LBGs 



The gravitational lensed LBG cB58 is the best-studied ty pical LBG. This ob ject is at z ~ 2.73 
and has a luminous n iass ~ lO^^Mp, a SFR ~ 40Moyr"^ (IPettini et all l2002f) and an effective 
radius of ri ~ 2kpc (|Seitz et al.lll998|) . for a 0.m = 0.3, = OJ,h = 0.70 cosmology. The 
average MIPS-LBG, with total mass ~ 10^ e stimated in this work, is more massive than cB58 



(cB58 is 10'" - 3 X IO'^Mq in iPipino et al.ll2011h . The model suggests a SFR ~ 200MQ/3;r for 



the average MIPS-LBG, which is higher than that of cB58. The metallicity of MIPS-LBGs (~ Zq) 
is higher than cB58 (~ 0.25Zo). These differences imply tha t MIPS-LBGs and cB58 are at the 



2010a|) and mass-metallicity (Maiolino et 



opposite ends of the observed mass-SFR (|Magdis et al. 
al., 2008) relations at z ~ 3. 

The shape of a galaxy SED depends on many physical processes as we have shown above. 
We compare the averag e MIPS-LBG in this work with other 24/i/Tj MIPS detected z ~ 3 LBGs 
(|Rigopoulou et al.ll2006l : IShim et al.ll2007h in Fig. [8]by plotting the estimated extinction E(V - B) 
versus observed 24/im flux density /24. As shown in Fig. [8l the observed frame does not 
correlate with the E{B - V) for these samples. One could not derive reliable IR properties only 
based on fitted parameters, which describe UV-optical properties. Therefore, we suggest that self- 
consistent SED models are needed when deriving detail properties of galaxies from their UV-radio 
multi-band data 



-20- 



3.6. Comparison with submillimetere galaxies (SMGs) 



The MIPS-LBGs with IR luminosity of L/« = 1.3 x IO'^Lq is in the class of ultra-luminous 
infrared galaxies (ULIRGs, Ljr > IO'^Lq). Since MIPS-LBGs could be the most massive, dusty, 
star-forming and young ellipticals as discussed in the above section, we suggest that these ellip- 
ticals contribute a significant population for ULIRGs at redshift z ~ 3. The non detection of 
infrared luminous LBGs in submm band is a puzzle, giv en their large and h igh UV-derived 



SFR, as well as substantial dust component (see detail in IChapman et alj (|2000|) and IShim et al 



feOOTh). There are only two detections of FIR count erparts of MIPS-LBGs reported in the litera- 
ture ([Chapman et aLlbOOOt IChapman & Caseyll2009|) . The higher dust temperature, different dust 
spatial distribu tion and lower SFR in LBGs than in SMGs have been s uggested to explain this 
challenge (e.g. IChapman et al.ll200d : iRigopoulou et alj|2006l : iMagdis et al.n2010b ). Based on our 
"best-fitting" model, our prediction of observational frame flux density of MIPS-LBGs at 850jum 
is 0.3mJy, which is under the confusion/detection limit of current submillimetere (submm) sur- 
veys (/sso = \mJy an d fs^n = 2m Jy). Note that, our 850/zm flux density /gso = Q3mJy is lower 
than /s5n = 1.36m7y ([Magdis et al.ll2010bh predicted by the "best-fitting" SED from templates 
(|Chary & Elbazll2001|) . while the IR luminosity of our model and other models are comparable 
(see Table 121). This reflects the fact that the predicted flux could be quite diff"erent at a particular 
wavelength, although a few observed fluxes or the overall-level SED (e.g. the IR luminosity) could 
be comparable, since the shape of SEDs could be quiet difl^erent. The detailed SED of a galaxy 
depends on all physical processes in our model. Besides the explanations mentioned in the liter- 
ature, we suggest the different galaxy morphology, dust environment and dust intrinsic properties 
could also cause the non detection of LBGs in submm blank surrey. However, it is hard to make a 
conclusion about the link between LBGs and SMGs before more data, such as a large population 
with detailed SEDs and chemical abundances, are available. 

The dust temperatures are merely those corresponding to the best-fitting single temperature 
grey-body model, and are not the actual dust temperature. We use the rest-frame lOO-to-850 //m 
flux density ratio versus 60-to-850 /im flux density ratio, /looZ/sso vs feo/fssa, as a proxy to FIR 
properties of galaxies. Our predictions of M310, Mil, M511 models from 0.1 Gyr to 0.7 Gyr 
are compared with the data of SCUBA local (recession velocity ~ few thousand km ) galaxies 
(|Dunne et al.ll2000|) in Fig. [9l Compared to SCUBA local galaxies, the young high redshift star- 
burst galaxies in our models have much higher /eoZ/sso given the same /looZ/sso- The more massive 
galaxies have higher /looZ/sso and /eoZ/sso- It is clear that the SCUBA local SMGs and young high 
redshift star-burst galaxies have different FIR properties. 



We argue that caution is needed when adopting the local empirical UV-IR relationships and 
templates for high redshift z ~ 3 LBGs to understand the physical processes in detail. 
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4. Conclusions 



In this paper, we modeled the rest frame UV-radio SED of the stacked MIPS-LBG at z ~ 3 
by a new "fitting" approach. In this self-consistent approach, we derived the average galactic-wide 
properties of MIPS-LBGs at z ~ 3 considering the stacked MIPS-LBG as a proxy of all MIPS- 
LBGs at z ~ 3. Our findings can be summarized as follows: 



1. The new "fitting" approach, which combines the chemical evolution model and GRASIL, 
can reproduce the rest frame UV-radio SED of the stacked MIPS-LBG at z ~ 3. This 
approach suggests that the MIPS-LBGs at z ~ 3 are likely young (0.3-0.6 Gyr) massive 
(~ IO^'Mq) elliptical galaxies with a fast star-burst regime of star formation. Chemical en- 
richment and dust evolution history are provided by chemical evolution models, calibrated 
at both low and high redshifts. 



We estimated that the the average stellar mass of MIPS-LBGs is in the range of ~ 6 x 10 



to ~ 1 X IO^^Mq. The "best-fitting" stellar mass and SFR of the stacked MIPS-LBG are 
8 X IO'^'Mq and 200A/o/jr, respectively. The stellar mass is in agreement with estimates 
done in other works and based on UV-optical SED fitting. This confirms that the derived 
stellar mass is robust in SED fitting. 

The "best-fitting" model of the stacked MIPS-LBG is the 0.5 Gyr lO^M© elliptical galaxy, 
which has dust mass of = 7 x 10^ and gas mass of A/^ = 1.6 x IO^^Mq. Our "best- 
fitting" dust mass of the stacked MIPS-LBG is about a factor of eight l ess than the value 



based on single temperature grey-body fitting (M^ = 5.5 + 1.6 x 10*^ M© in lRigopoulou et al 



2010h . 



The parameters of molecular clouds and dust of the Milky Way cannot fit the stacked MIPS- 
LBG SED. Our "best-fitting" parameters reflect that: i) MCs of the stacked MIPS-LBG are 
more dense dusty environments than in the Milky Way; ii) the dust size distributions in 
stacked MIPS-LBG may be flatter than in the Milky Way; iii) both small and big size dust 
are needed to reproduce the stacked MIPS-LBG SED; iv) non-negligible PAH also make an 
important contribution to MIR flux. More observational data, such as chemical abundances 
and emission lines, are needed to make a more realistic prediction of the properties of all 
MIPS-LBGs, especially for dust intrinsic properties. 

We suggest that high redshifts star-forming ellipticals make a significant contribution to 
the population of ULIRGs. The dust properties and morphology of high redshift star-burst 
galaxies could be different from local star-burst galaxies. We argue that self-consistent SED 
models are needed to investigate the detailed properties of high redshift LBGs . 
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Fig. 5. — Rest-frame SEDs of a lO'^M© galaxy at 0.5 Gyr with various dust environments and 
dust intrinsic properties. Blue lines: various dust size distributions and MC properties. Only with 
these various properties, the UV-optical part of the SED varies by a factor of ten. It is clear that the 
MW properties cannot fit t he SED. Red l ine: MW dust and MC properties in Tableland Table |4] 
MIPS-LBGs data are from lMagdis et all (l2010bh as described in FigS 
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Fig. 6.— "Best-fitting" rest-frame SED. Red solid line: the "best-fitting" SED of a lO^M© galaxy 
at 0.5 Gyr with MIPS-LBGs dust and MC properties shown in Tableland Table IH At the third 
fitting step, this "best-fitting" model is given by min-;^^^ fitting approach based on the grid models 
shown in Fig. [51 The values in Table [3] and Table H] just give the general properties of the dust and 
MCs and should not be considered as the real ones. Note that, one can get an another "best-fitting" 
model using another galaxy model (such as a IO^^Mq galaxy at 0.4 Gyr). However, the overall- 
level SEDs are similar for models within a small age parameter space (see Fig. H]). Therefore, the 
trend of the dust and MC properties is the same for those models with different age. However, 
the parameters of Milky Way can not fit the SED no matter what the galactic age is adopted (see 
examples in Fig. S]). The mean differences between the "best-fitting" and "Milky Way parameters- 
fitting" SED (Blue dashed line, "W MP-fitting") are the UV -optical shape and the position of the 
FIR peak. MIPS-LBGs data are from lMagdis et al.l JlOlObh as described in FigS 
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Fig. 7. — Rest-frame SEDs of different models. To show the necessary of more reliable SFHs, 
SEDs of a model (age ~ 0.5Gyr, stellar mass ~ 8 x 10'*' M©), which takes a self-consistent SFH ( 
Mil case in Fig. [U SSF for short) and a model (age ~ 2G yr, stellar mas s ~ 9 x lO'^M©) with a 



constant SFH (CSF) are shown. MIPS-LBGs data are from lMagdis et al.l (l2010b|) as described in 
Fig [31 
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Fig. 8. — Observed 24 /jm flux density vs estimated extinction E(B-V) for 24/zm M I PS de - 
tected LBGs at z ~ 3. Asterisk, cross, square symbols represnet data from IShim et al.1 (|2007|) . 



Rigopoulou et al.l (|2006h and the average MIPS-LBG, respectively. The f24 is not correlated with 
E(B-V) for these objects. 



-31 - 




200 400 600 800 1000 



60 850 



Fig. 9. — Rest-frame finn/ fssn flux density ratio vs /eoZ/sso flux density ratio. The data of SCUBA 
local galaxies are from lDunne et al.l(|2000|) . Lines are our model predictions of /looZ/sso vs /eoZ/sso 
evolution from 0.1 to 0.7 Gyr. The large triangle, asterisk and square are corresponding for M310, 
Mil and M510 model predictions at 0.5Gyr, respectively. The high redshift galaxies and local 
SMGs have different FIR properties. 



